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ABSTRACT 

We present the 24 /zm rest-frame luminosity function (LF) of star-forming galaxies in the redshift 
range 0.0 < z < 0.6 constructed from 4047 spectroscopic redshifts from the AGN and Galaxy Evo- 
lution Survey of 24 /im-selected sources in the Bootes field of the NOAO Deep Wide-Field Survey. 
This sample provides the best available combination of large area (9 deg^), depth, and statistically 
complete spectroscopic observations, allowing us to probe the evolution of the 24 /im LF of galaxies 
at low and intermediate redshifts while minimizing the effects of cosmic variance. In order to use the 
observed 24 /im luminosity as a tracer for star formation, active galactic nuclei (AGNs) that could 
contribute significantly at 24 /xm are identified and excluded from our star-forming galaxy sample 
based on their mid-IR spectral energy distributions or the detection of X-ray emission. Optical emis- 
sion line diagnostics are considered for AGN identification, but we find that 24 /im emission from 
optically-selected AGN is usually from star-forming activity and therefore should not be excluded. 
The evolution of the 24 /im LF of star-forming galaxies for redshifts of z < 0.65 is consistent with 
a pure luminosity evolution where the characteristic 24 /im luminosity evolves as (1 -I- 0)3-8±O-3_ -y^g 
extend our evolutionary study to encompass 0.0 < z < 1.2 by combining our data with that of the 
Far-Infrared Deep Extragalactic Legacy Survey. Over this entire redshift range the evolution of the 
characteristic 24 /tm luminosity is described by a slightly shallower power law of (1 -I- z)^-'^^'^-^. We 
find a local star formation rate density of (1.09 ± 0.21) x 10~^ Af0yr~^Mpc^'^, and that it evolves as 
(1 + z)^-^^^-^ over 0.0 < z < 1.2. These estimates are in good agreement with the rates using optical 
and UV fluxes corrected for the effects of intrinsic extinction in the observed sources. This agreement 
confirms that star formation at z < 1.2 is robustly traced by 24 /im observations and that it largely 
occurs in obscured regions of galaxies. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies 



1. INTRODUCTION 

We have learned in the past 15 years that the star 
formation rate (SFR) of the universe peaked at a red- 
shift of z > 1 and then rapidly declined by an or- 
der of magnitude to the present day (e.g.. ILillv et al.| 
1996; Ma dau et all '1996: Chary & Elbaz 2001; Hopkindl 
2004: Le Floc'h et al, 2005 : Perez-Gonzalez et al. ,2005i; 
Hopkins fc BeacomI |2006[) . The exact cause of the de- 
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cline of the SFR remains unknown. However, since the 
period from z ~ 0.6 to the present encompasses nearly 
half of cosmic time coinciding with the SFR decline, this 
redshift range is a key regime for constraining the mech- 
anisms that may be responsible. 

The evolution of the SFR has been measured exten- 
sively at mid-IR (MIR; 8-30 /im), far-IR (FIR), ultravi- 
olet (UV), and optical wavelengths. Among these wave- 
lengths, MIR and UV are ideal tracers for dust-obscured 
and unobscured galaxies, respectively. The MIR emis- 
sion of galaxies originates from the reprocessed UV ra- 
diation of dust-embedded, highly obscured star-forming 
regions, and hence is sensitive to radiation from star for- 
mation that can not be seen directly at optical wave- 
lengths (e.g., Eiinicutl [TpSl; |Buai^^^ The 
stellar UV radiation field is largely created by hot OB 
stars (M* > 10 Mq), which have very short lifetimes 
(< 20 Myr) and hence trace current star formation. The 
dust obscuration is generally substantial for star forma- 
tion in galaxies at and above the characteristic luminosity 
L* (Schmit t et al..,2006.: Buat et al . 2007). Such galaxies 
comprise the majority of galaxies studied beyond the low 
redshift volume, so the MIR offers a relatively complete 
probe of SFR at intermediate and high redshifts. Among 
the IR wavelengths, the 24 /im-band is particularly suit- 
able to study SFR evolution because it avoids most of the 
strong emission features from the stochastically heated 
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aromatic hydrocarbons (hereafter PAH) at wavelengths 
3.3 < A < 13 /Ltm for redshifts z < 0.5 while still avoiding 
the contamin ation from the Ga lactic IR cirrus at longer 
wavelengths (jShupe et al.lll998[ ). 

Measurements of IR number counts and luminos- 
ity functions (LFs) of galaxies demonstrate the rapid 
evolution of the SFR. Early IR studies utilized the 
IRAS Bright Galaxy Sample (iSo ifer et all 119891) and 
the IRAS calibration field (Hac king fc Houckl fl987l) 
and focused on accounting for the cosmic infrared 
background (GIB; Elbaz pt all 119991) with di screet 
MIR sources (IHackiiig et al.lll987t iSaunders et al.lll990l : 
ISoifer fc Neugebaueii Il99ii r Studies of galaxy num- 
ber counts with ISO and Spitzer have shown that 
IR galaxi es evolv e very strongly ('Elbaz e t al.l 119991: 
Franc cschi ni et al.l[20 01:'D ole et all [2004: Ma r lean et ail 
20041: iPapovic h et al.ll2004l ). As redshifts of IR samples 
became availabl e, th e LF could be constructed i n the 
nearby volume ([ Riekc & Lcbofskv 1986; Saunders et al.l 
19901: lYahil et alj 11991; .Rowan- Robinson et al.. 11997 



Shupe et al.l Il998l i and subsequently extended to inter- 
mediate and high reds hifts (z > 3) with large samples of 
pho tometric redsh ifts (iPozzi et al . 2004; Le Floc'h et ah 



2005; Pcrcz-Gonzalez et al 



WmA iGaputi et all 120071 : 



Marleau et al..,.2007i ) to study the SFR evolution over cos- 
mic time. However, most of these studies were either in 
the local regime (z < 0.1) or at intermediate to high 
redshifts (z > 0.5). The redshift range 0.1 < z < 0.6 
includes more than 4 Gyr of galaxy evolution and needs 
to be characterized in detail to understand this process. 

Observations of galaxy evolution at low and interme- 
diate redshifts face two challenges. First, a large area, 
deep survey is required to access a large volume at lower 
redshifts, which is critical to minimize cosmic variance. 
For example, a deep survey at redshift z ~ 0.5 would 
need five times the solid angle coverage of one at z ~ 1 
to probe a comparable comoving volume. Second, while 
photo-z's are generally adequate for distance determina- 
tion at high redshift, accurate spectroscopic redshifts are 
critical at low and intermediate redshifts. To illustrate 
this issue, consider typical photometric redshift uncer- 
tainties of Az/(1 + z) — 0.05. The corresponding frac- 
tional uncertainty in z would be a tolerable 10% at z 1, 
while at z ~ 0.25 and below, the fractional uncertainty 
would exceed 25%. Also, photometric redshifts are only 
accurate for galaxies whose spectral energy distributions 
(SEDs) are well matched by either SED templates or, in 
the case of empirical photometric redshifts, galaxies in a 
spectroscopic training set. This may not be the case for 
heavily obscured galaxies such as luminous IR galaxies 
(LIRGs). In addition, the accuracy of some photometric 
redshift samples in the literature has not been verified 
with spectroscopy, and unknown systematic errors may 
be present. 

In this work we construct the rest-frame 24 /im LFs of 
star-forming galaxies, measure their evolution, and esti- 
mate the SFR for 0.0 < z < 0.6 from the LF. In addition 
to constraining the LF and SFR evolution accurately at 
low and intermediate redshifts, we combine our results 
with those from the Far-IR Deep Ex tragalactic Legacy 
Survey (FIDEL, iMagnelli et"aIll2009D to extend the con- 
straints to z ~ 1.2. We organize this paper as follows. 
In Section [2] we discuss our Spitzer 24 fim sample, the 
spectroscopic redshifts from AGN and Galaxies Evolu- 



tion Survey (AGES), the optical catalog from NDWFS, 
and address their completeness. In Section|3]we describe 
the construction of the LF, starting with the derivations 
of galaxy luminosities (Section l3.ip , identification and ex- 
clusion of AGN from star- forming galaxies (Section (32]), 
estimators for the LF (Section 13. 3|) . parameterization of 
the evolution of the LF with redshift (Section 13. 5p . and 
the estimation of the SFR from the LF (Section 13. 6p . 
We then present the results in Section |4] for the local LF 
(LLF; Section |4T|) . the LF evolution at low and interme- 
diate redshifts (Scction l4.2|) . the LF evolution constraint 
out to z '--^ 1.2 (Section l4.3p . as well as the SFR evolution 
of star- forming galaxies (Section 14. 4p . We summarize 
our results in Section [5j We assume a ACDM cosmology 
with = 0.3, 17a = 0.7, and Hq = 70 km s-^Mpc'^ 
throughout this paper. 

2. THE DATA 

We construct the 24 /im LF of galaxies in the 9 sq. de- 
gree Bootes field of the NOAO Deep Wide-Field Survey 
(NDWF S Data Release 3; a = 14''32™, S = -h34°16', 
J2000.0; lJannuzi fc Devlll999[ ). The Bootes field has op- 
tical imaging reaching the 5-a flux limit of 26.6, 26.0, 
26.0, and 21.4 in Bw, R, I, and K (Vega magnitudes), 
respectively. Details a bout the NDWFS can be found in 
iJannuzi fc Devi (|1999D and at the survey 's Web siteQ. 

In addition to the deep optical imaging, the Bootes 
field has extensive panchromatic coverage from other sur- 
veys, including an X -ray survey with the Chandra Space 
Telescope (XBo6tes: lMurrav et al.ll2005D comprising of 5 
ks AGIS- I exposures, near-U V deep 30 ks imaging from 
GALEX (|Martin et al.ll200l . the IRAG Shallow Survey 
in the 3.6, 4.5, 5.8, and 8.0 fj,m bands reaching mag- 
nitudes of 19 .1, 18.3, 15.9, and 15. 2 mag (Vega mag), 
respectively, (jEisenhardt et al.ll2004D the deeper Spitzer 
Deep, Wide-Field Survey (SDWFS. lAshbv et al.l 12009') . 
and MIPS MIR photometry at 24 down to 0.27 mJy 
(jPapovich et al.l 120041; iHouck et al. 2005). These imag- 
ing data are complemented by the set of optical spec- 
troscopic redshifts from the AGN and Galaxy Evolution 
Survey (AGES; Kochanek et al., in prep.), which is a sta- 
tistically complete redshift survey of the field at low and 
intermediate redshifts (0 < z < 0.8) (C. S. Kochanek et 
al., in preparation). 

2.1. Spitzer 24 ^m- Observations, Reduction, and the 
Final Catalogue 

The 24 /im photometry in the Bootes field was obtained 
by the Spitzer Space Telescope wide-area survey with 
the Multiband Im aging Photometer for Spitzer (MIPS; 
IRieke et all l200l and proc essed with the MIP S GTO 
Data Analysis Tool (DAT; iGordon et all 120051) . The 
mean exposure time per pixel across the Bootes field is 
87 s. The pixel size of the final mosaic is 1'.'25 with a 
FWHM of 5'.'7. 

Source extraction from the final 24 / im mosaic was 
done using the techniques described in IPapovich et al.l 
(I2OOI . Briefly, we weight the image by the typical inte- 
gration time of each region to obtain uniform noise across 
the mosaicked image. We filtered the image with a Gaus- 
sian kernel similar to our 24 /im point-spread function 

^ http:/ /www. noao.edu/noao/noaodeep 
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(PSF) to remove artifacts and opti mize point-sou rce de- 
tection and then used DAOPHOT ()Stetsonlll987l ) to de- 
tect sources. The actual photometric measurement was 
carried out by fitting empirical PSF constructed from 
bright sour ces in the image . Base d on Monte Carlo sim- 
ulations bv lPapovich et all ()2004D the 80% completeness 
flux threshold is F,y(2A /im) = 0.27 mJy, and there are 
8,572 objects in our sample above this flux level. 

The MIPS survey of the Bootes field combines con- 
siderable depth with a large contiguous area of 32,457 
arcmin^ (9.0 deg^). While many other fields have 
been surveyed to deeper F^{2A /x m), such as the 
GOODS, ELAI S and CDF-S fields (IPozzi et al l [2001 
ILe Floc'h etan [200l ICharvl [20071: lElbaz et all [2007f) 
with 80% completeness flux of 25 fi Jy, 61 fj, Jy and 83 /iJy, 
respectively, they also cover significantly smaller area 
(0.011 deg2, 0.036 deg^, and 0.407 deg^, respectively), 
and therefore encompass significantly smaller volumes 
at low and intermediat e redshifts. The Co smic Evolu- 
tion Survey (COSMOS, IScoville et all [2001 h as a large 
sample of spectroscopic redshifts (iLillv et al.ll2007D an d 
deep, -80 fiJy, 24 ^m imaging (jLeFloc'h et al.l 120091 ). 
The 2 def area of the COSMOS field is one of the best 
combinations of large area and depth at higher redshifts; 
however, a larger angular coverage is still required to 
suppress the effects of cosmic variances at low and inter- 
mediate redshifts. The available large-area MIR surveys, 
such as the 5piteer Wide-area Infrared Extragalactic Sur- 
vey (SWIRE) and the European Large Area ISO Survey 
(ELAIS, 15 /im), currently lack thorough spectroscopic 
redshifts. 

2.2. Optical Spectroscopy Target Selection, Observation, 
and the Final Catalogue 

The Bootes field was observed spectroscopically by the 
AGES during 2004-2006 using the Hectospec spectro- 
graph in a survey region of 7.9 deg^ (C. S. Kochanek et 
al., in preparation). Hectospec is a fiber-fed spectrograph 
with 300 robotically assigned fibers covering a 1-degree 
diameter field of vie w at the //5 focus of the 6.5-meter 
MMT Observatory (|Fabricant et al.ll2005D . AGES pro- 
duced -20,000 redshifts, of which -16,500 and -3,500 
are galaxies and quasars, respectively, with redshift un- 
certainties less_than_a^_— 100 km (as in the similar 
survey bv iPapovich et all ([2006)) and a median galaxy 
redshift of z = 0.26 (C. S. Kochanek et al., in prepara- 
tion). 

Targets for spectroscopy were principally selected on 
the basis of their optical magnitudes. AGES targeted 
galaxies to /c < 20, where Ic is the Kron-likc magnitude 
(in Vega units) in the Kron-Cousins I-band of the ND - 
WFS (SExtractor's MAG_AUTO: [Bertin fc Arnoutslil99i ). 
For galaxies brighter than Ic = 18.5, there was full sam- 
pling. For fainter galaxies with 18.5 < Ic < 20, AGES 
used random sparse sampling of 1 in 5 galaxies, with 
rates that depended on a variety of other factors, so as 
to put emphasis on more unusual objects. Additional 
targets were added to the spectroscopic sample if they 
were brighter than certain fiux limits in the X-ray, in- 
frared, other optical bands, and radio. For example, all 
galaxies with F,y(24 //m) above 0.5 mJy were targeted, as 
well as 30% of those down to 0.3 mJy. For each object, 
the a priori sampling rate is known and so the object 



can be exactly weighted to form an unbiased sample of 
Ic < 20 galaxies. AGES observed the Bootes field with 
15 complete tilings, so the completeness for the target 
galaxies is very high. Our final sample has 4.3% incom- 
pleteness due to difficulties in assigning fibers to targets 
and 2.1% incompleteness due to failing to get a success- 
ful redshift from an assigned fiber. We also estimate that 
the parent photometry catalog is missing 4% of the ob- 
jects in the range 15 < Ic < 20, half of which is due to 
interference by nearby bright stars, the rest likely due to 
blending and corruption of the parent photometry. We 
model these mild incompletenesses, including factors for 
areal target density and surface brightness, and correct 
the final weights of objects in the redshift catalog (R. J. 
Cool et al., in prep.). As the total incompleteness cor- 
rection is only 10%, it is very likely that the systematic 
error in the correction is far less than the Poisson and 
cosmic variance in our results. 

Each Hectospec fiber is 1'.'5 in diameter and produces 
a spectrum covering the wavelength range of 3500 A 
< A < 9200 A with 6 A FWHM resolution {R - 1,000). 
Typical signal-to-noise ratios of the AGES spectra are 
2 — 14 per pixel with a median of 5 per pixel. We fit 
redshifts to AGES with two automated pipelines, one 
based on the SDSS template-fitting code, the other on 
cross-correlations. The pipelines make use of spectral 
features across the spectral range (absorption lines, emis- 
sion lines, and spectral breaks due to metal-line ab- 
sorptions) and hence result in accurate redshifts regard- 
less of whether a particular spectral line is compromised 
(e.g. by overlap with sky lines). The two pipelines gave 
broadly consistent results. In addition, we visually ex- 
amined every fit. Any objects with visually questionable 
results were tagged for reobservation, since AGES cov- 
ered the full field in six sequential passes. The result is 
a spectroscopic catalog with high completeness and high 
reliability. The spectroscopic success rate, after fiber as- 
signment, is 98%, and we expect the redshift error rate on 
claimed successes to be well below 1%. Any systematic 
biases due to errors in the completeness corrections or 
the redshift catalog are likely to be considerably smaller 
than the statistical errors of our LF. 

In the redshift range of interest, 0.05 < z < 0.65, we 
have redshifts for 4,047 objects with Fi,{24 fim) above 
the 80% completeness limit of 0.27 mJy. 

3. INFRARED LUMINOSITY FUNCTION 

Here we describe the construction of the 24 /im LF. 
We start by estimating the 24 /im and optical (I- 
band) rest-frame luminosities, then identify and exclude 
AGN in our sample, and finally construct the LF using 
the c lassical 1/Knax ([Schmidt! 119681 iHuchra fc Sargentl 
| 1973fl and the m aximum likelihood estimator (MLE; e.g., 
'Sand age et all 119791 [Marshall etlll [1983[ : lYahil et alJ 
1991). Later in this section, we describe the parame- 
terization of LF evolution, which will be used to study 
the evolution of the SFR. 

3.1. 24 jJ-m and I-band Luminosity Estimation 

In addition to the obvious need for 24 /im luminosi- 
ties, L(24 /im), to construct the LF, we have to com- 
pute the I-band luminosity for each galaxy. AGES spec- 
troscopy was targeted based on an I-band, flux-limited 
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parent catalog. Hence, an object can only be included in 
the final LF sample if it is detected above the flux limits 
at both 24 /im and I-band. Therefore, the luminosities 
in both bands are needed for estimating the maximum 
redshift, Zmax, that an object would still be detected in 
both bands to the construct the 1/V max LF. L(24 /xm) i s 
estimated by the formalism given bv lRiekeet"£a1(|200a ). 
For the I-band, we derive the I-band k-correc tion using 
the kcorrect code (jBlanton fc RoweisI |2007() and then 
use the k-correction to estimate the I-band luminosity. 

The I-band k-corrections were derived using optical 
photometry from NDWFS and redshifts from AGES. A 
fraction of our I-band objects are affected by scattered 
light from bright stars. Specifically, some of the NDWFS 
I-band Kron-like magnitudes (SExtractor's MAG_AUTD), 
Ic, are overly bright due to contamination from low- 
surface-brightness halos of nearby bright stars. This 
problem can be corrected by synthesizing a new I-band 
magnitude, Ir, by adding a Kron-like magnitude in R- 
band. Re, to the I — R color from 6" apertures {In = 
Rc + Ig" — R&' )- This synthesized Ir is less sensitive to 
the I-band halo than Ic- Then we compare Ir to Ic to 
compute a new I-band magnitude, Itotai, that uses the 
fainter magnitude if the two differ significantly; otherwise 
we average them (R. J. Cool et al., in preparation). The 
hotai computation is done in a way that the continuity 
of the magnitudes is preserved. In the cases that require 
averaging Ir and Ic, the differences are on the order of 
a tenth of a magnitude and hence averaging a logarithm 
does not introduce a significant bias. We note that the 
aperture correction is insensitive to the choice of estima- 
tion procedure and that the final LFs are little affected by 
this correction, ultimately because the I-band generally 
has deeper photometry than at 24 /im and hence most of 
the values of z„iax will be limited by the 24 /im photom- 
etry. Apart from the Itotai for I-band photometry, we 
use the 4"-diameter aperture magnitude from NDWFS 
for the Bw and R bands to perform k-correction esti- 
mation. We adopt the following corrections to convert 
the Vega zero points to AB: +0.02 mag, -1-0.19 mag, and 
-1-0.44 mag in Bw, R, and /, respectively. The Galactic 
extinction for our field is Aj = 0.02 mag ([Schlegel et al.l 
[19981) . 

The estimation of L(24 /im) requires an assump- 
tion for the spectral energy distribution (SED) of a 
galaxy in the IR reg ime, which has been modeled by 
many authors fe.g. iDale et al.l 120011: iCharv fc Elba: 
2001; Dale & Hclo ul 120021: iLaeache et al.l 120031 1200 
&mith et al. ,2007. 
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Dale et al.l 120071) 



More recently, 
constructed a library of SEDs by 
combining stellar photospheric spectra from 0.4—5 /tm 
with Spitzer IRS spectra from 5—30 /tm and a mod- 
ified blackbody spectrum at longer wavelengths for 
11 representative local Luminous IR galaxies (LIRGs; 
10"io < L(TIR) < IO^^Lq) and ultraluminous IR 
galaxies (ULIRGs; L(TIR) > IO^^Lq). Higher L(TIR) 
galaxies contain more intense star-forming activity and 
hence a larger fraction of warm dust, leading to changes 
in the SED with IR luminosity. 

Normally, the derivation of SFR from MIR photom- 
etry involves selecting the SED that best matches the 
object's luminosity given its MIR flux and redshift, inte- 
grating this SED to obtain i(TIR), and then using a SFR 



metric to calculate the SFR. However, this procedure im- 
plies fitting a SED with only one data point, 24 /tm in 
our case, because the number of objects that are detected 
in multiple MIPS bands is ve ry limited, which c ould lead 
to significant uncertainties. IRieke et al.l (|2009[ ) circum- 
vented this problem by first assigning a SFR to each 
SED, which has a unique L(TIR). They then derive a 
flux in the relevant band pass (24 /tm) over a range of 
redshifts. 

Giv en that i(TIR)/L(24 /im) ~ 10 (e.g.. IRieke eFal] 
I2009D . our sample with galaxies mainly in the 24 /im lu- 
minosity range 10^-^ < i(24 /im) < 10-^°'^ (Figure [J), is 
dominated by normal star-forming galaxies (L(TIR) < 
lO^^L©). Figure [1] also illustrates that LIRGs only start 
contributing significantly to the MIR luminosity density 
at redshifts greater than z = 0.2. Although ULIRGs 
start to appear in our sample at z ^ 0.4, we will demon- 
strate in Section 13.2.41 that they are mostly AGN and 
we will exclude them from our sample of star-forming 
galaxies. 

3.2. Separating Star-Forming Galaxies and AGN Hosts 

To quantify the 24 /tm emission due to star formation, 
we need to minimize the contribution of AGN. There are 
several methods to identify AGN in a galaxy sample. Op- 
tical emission l i ne di agnostics such as the BPT method 
([Baldwin et al.l [l98l[ ) are commonly used when optical 
spectra are a vailable. MIR power -law characte ristics 
of AG N Ce.g.. lAlonso-Herrero et al.l 12006: Donley et al.l 
20071) and MIR color criteria (e.g., iLacv et al.i 12004 



Stern et al.l [20051) can be used if multi-band MIR pho- 
tometry is available. X-ray emission is also a direct in- 
dicator of AGN activity. Since the Bootes field has mul- 
tiwavelength coverage that permits using all these, we 
investigate all of them in this section. 

3.2.1. Optical Emission Line AGN Diagnostics 

In addition to identifying the broad-line AGN (type 1), 
our spectroscopic survey allows us to identify narrow- 
line AGN (type 2) in the sample using the BPT dia- 
gram. We require that all four relevant lines (Ha, H/3, 
[OIII], and [Nil]) are detected at S/N > 3. Thus, the 
BPT method can only be used for redshift of z < 0.35 
because at higher redshifts the Ha emission line shifts 
beyond the red limit of the spectral range covered by 
Hectospec. We identify AGN on the BPT diagram (Fig- 
ure [21 left) as objects with line ratio s above the theoret- 
ical h mit (hereafter the Kewley line; iKewlev et al.|[200ll 
I2006D of ratios that star-forming regions could produce. 
The Kewley line was estimated by considering all pos- 
sible combinations of star-forming galaxy spectra based 
on the Starburst99 (jLeitherer et al.l ll999) models. Ob- 
jects above the Kewley line will be referred to hereafter 
as optically selected AGN. We identified 421 optically se- 
lected AGN at 0.0 < z < 0.35, 46% of which (192/421) 
have Fj.(24 /im) > 0.27 mJy. 

However, we need to determine whether the 24 /im 
emission from optically selected AGN is actually domi- 
nated by the AGN. Our 1'.'5 diameter spectroscopic fiber 
would cover a distance of 6 kpc at the median redshift 
of the AGES sample {z = 0.26) and the fiber would be 
placed at the center of the target galaxy. Thus our BPT 
classification would be sensitive to the nuclear flux from 
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Figure 1. Redshift distribution of the rest-frame 24 /im luminosity for our sample. The gray dots show all 24 fim galaxies with spectroscopic 
redshifts from AGES. The black dots are those above the 24 fim completeness flux limit (0.27 mjy) at 0.05 < z < 0.65 that we use to 
construct the LF. The line indicates the L(24 fim) corresponding to this completeness flux threshold as a function of redshift. Inset: the 
redshift distribution of the AGES objects detected at 24 fj,m. 
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Figure 2. Left: optical selection of AGN using the BPT diagram. Orange diamonds and magenta dots are power-law and X-ray selected 
AGNs, respectively. The Kewley and Kauffmann lines are shown by red dot-dashed line and blue dashed line, respectively IIKewlev et al.l 
120011 : rRauffmann et al. 2003). Many power-law and X-ray selected AGN are identified as star-forming galaxies by the BPT diagram. Note 
that a large fraction of power-law and X-ray selected AGN are not shown here because they do not have sufficiently strong emission lines 
to apply the BPT diagnostic and most of them are at z > 0.35 where Ha line redshifts out of the Hectospec spectrum. Right: the Venn 
diagram for 712 AGNs at 0.0 < z < 0.65 to illustrate the overlapping of AGN identified by X-ray detections, power-law fitting, and the 
UStern et al.l 2005) color criteria. The numbers in each area and the numbers in parentheses are the total numbers of objects and the 
numbers of objects detected at 24 fim, respectively. 288 of 712 AGNs in this Venn diagram are detected at 24 fim, pass all quality cuts, 
and are in the main galaxy sample. They are excluded from the sample to obtain the star-forming galaxy sample that can be accurately 
used as a tracer of star formation. 



the AGN but could miss the emission hues contributions 
from star-forming regions outside the fiber. In the case 
that star formation dominates the global optical emission 
but an AGN dominates the emission within the fiber, 
the global 24 /xm flux could in fact be dominated by the 
emission from star-forming activity and excluding such 
objects would lead to an underestimation of the star for- 
mation rate. 

We therefore investigate the origin of mid-IR emissions 
from the optically selected AGN by studying the IRAC 
colors of these AGN. Using IRAC photometry from the 
IRAC Shallow Survey (6" aperture, Section[2|), we study 
the [3.6 fim\ — [4.5 fim\ versus [5.8 fim\ — [8.0 fim\ color- 
color diagram (hereafter, [3.6] — [4.5] and [5.8] — [8.0], re- 
spectively) for optically selected AGN (Figure El left). 
The color-color diagram can be used to identify MIR 
aromat ic emission chara c teristic of star formation ac- 
tivity (|Smith et al.l l2007t iShi et al.ll2007[ ). because the 



strongest aromatic features at 6.2, 7.7, and 8.6 /im fall 
into the 8.0 /im band (channel 4; ~ 6.2 — 10 /im) at 
redshifts 0.0 < z < 0.5. This behavior results in a redder 
[5.8] — [8.0] color for objects with MIR aromatic emission, 
which also implies that the MIR emission of these ob- 
jects is dominated by star formation. The empirical color 
boundary that separates the arom atic and non-aromatic 
region is [5.8] -[8.0] ~ 1 (|Brand et al. 2009). We found 
optically selected ACN distributed across the full range 
of [5.8] — [8.0] color. However, if we further separate op- 
tically selected AGN into those detected at 24 /zm and 
those that are not, they split into two locales. Optically 
selected AGN without 24 /im detection are relatively blue 
in [5.8] — [8.0] color while more than 70% of those with 
24 /im detection congregate in the aromatic locus. The 
bluer MIR locus of optically selected AGN without 24 /im 
detection is consistent with the Rayleigh- Jeans tail of an 
old stellar population's photospheric emission that peaks 
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Figure 3. Mid-IR, properties of AGNs in our sample. Left: tiie [3.6] — [4.5] vs. [5.8] — [8.0] color-col or diagram separates objects with and 
without aromatic emissions at an empirical boundary of [5.8] — [8.0] ^ 1 (vertical line. lBrand et ani2009f) . with aromatic-emitting objects 
on the redder [5.8] — [8.0] side. Blue and red dots represent optically selected AGN with and without 24 fim detection, respectively. Right: 
the same diagram showing the IStern et al.l (2005) AGN selection criteria (trapezoidal box), power-law selected AGN (orange diamonds), 
and X-ray selected AGN (magenta dots). A majority of optically selected AGN that emit strongly at 24 fim are aromatic-emitting objects 
and their MIR emission is likely dominated by star-forming activity. On the other hand, a majority of X-ray selected AGN share the same 
locus as power-law selected AGN, indicating that their AGN component could contribute significantly at 24 fim. 



at 1.6 /im. This suggests the source of 24 /im emission 
in an opticahy selected AGN is hkely to be star-forming 
activity, not the AGN. 

3.2.2. X-ray Detections 

We have also searched for AGN detect ed in the 5 ks 
X-ray survey of the Bootes field (XBootes; lIVTurrav et al.l 
[20051) . We define an X-ray selected AGN as any source 
with two or more X-ray counts and a > 25% probabil- 
ity based on t he Ba yesian matching method described by 
iBrand et al.l ([200l . Since XBootes is designed to have 
a large contiguous area and a shallow flux limit, X-ray 
identification based on the survey will only identify the 
strongest AGN. The X-ray luminosities of these sources 
are t ypically brighter than IC*^ ergs~^ ([Hickox et al.l 
l2009f ) and hence they are unlikely to be dominated by 
star-forming activity. Although it is not possible to infer 
directly the AGN contribution to the 24 fim flux from the 
strength of the X-ray emission, stro ng AGN are known t o 
emit significantly in the IVIIR (e.g., Barmbv et alj [20061 ) . 
We found 498 X-ray selected AGN at 0.0 < z < 0.65 and 
175 of them are at the redshift range where we construct 
the LF (0.05 < z < 0.65), are in the main galaxy sample, 
and are detected at 24 ^m. 

3.2.3. MIR AGN Diagnostics 

AGN with significant IVIIR emission can be identified 
most directly by applying IVIIR selection methods. IVIIR 
photometry can be used to identify luminous AGN due 
to their characteristic power-law SED rising from the 
optical wavelengths to the IR. This power-law SED is 
due to the combination of thermal emission of warm 
and hot dust comp onents and non-thermal emission of 
the nuclear region CNeugeb auer et al.l [TOTOt lElvis et al.l 
[1994; Riekc & Lcbofskv 19811). AGN can be identified 
by fitting a power law, = i/", to all four bands of 
IRAC photometry covering 3.6-8 fj,m (jEisenhardt et all 
I2004D with minimization. We require an object to 
be detected in all four bands at above 5-cr for power- 
law fitting. Power-law selected AGN are defined as 



galaxies e xhibiting power-laws with a ranging from —0.5 
to -2.0 dlv ezic et al.l [20021: lAlonso-Herrero et al.l [20061: 
IDonlev et a l. 2007). The power- law method identified 
169 AGN. Among these AGN, 144 are detected at 24 fim 
(85%), which is the highest fraction among the AGN clas- 
sification methods. This suggests that the 24 fj,m emis- 
sion of power-law selected AGN is an extension of their 
power-law AGN SED to longer wavelengths rather than 
from the warm dust characteristic of star-forming galax- 
ies. Therefore the power-law selection of AGN provides a 
good basis to exclude AGN from our 24 /j,m star-forming 
galaxy sample. However, the signal-to-noise requirement 
on the IRAC photometry for power-law fitting affects the 
completeness of power-law AGN selection. Among these 
169 AGN, 73 are at 0.05 < z < 0.65 and pass all target- 
ing and quality cuts, and are detected at 24 /im. 

Another manifestation of power-law SED of AGN is 
their locus on IVIIR color-color diagrams (Section I3.2.ip . 
Again, we consider the [3.6] — [4.5] vs. [5.8] — [8.0] color- 
color diagram. AGN (mostly type 1) exhibit power-law 
SEDs and this results in a distinct locus that is signifi- 
cantly redder in the [3.6] — [4.5] color than bot h early type 
and star-forming galaxies (see Figure 1 of IStern et all 
|2005[) . The locus of early type galaxies, which also ex- 
hibit power-law SEDs but with a > —0.5, would lie 
blueward in both color indices (lower-left locus on the 
diagram) because, as mentioned earlier, their SEDs in 
the NIR are domin ated by the Rayleig h- Jeans tail of old 
stellar populations. IStern et al.l ([2005^ empirically define 
the AGN locus on the color-color diagram (the "Stern 
wedge". Figure [11 right). AGN selected based on this 
color-color criteria (hereafter IRAG color selected AGN) 
are defined by galaxies with ([5.8] — [8.0]) > 0.6 and 
([3.6] -[4.5]) > 0.2([5.8]-[8.0])-K0.18and ([3.6] - [4.5]) > 
2.5([5.8] - [8.0]) - 3.5. We found 354 IRAC color se- 
lected AGN in the redshift range 0.0 < z < 0.65, 80% 
(282/354) of which have 24 /im detections. A major- 
ity (86%, 146/169) of the powe r -law selected AGN are 
also selected by the IStern et al.l (120051 ) criteria. The re- 
verse is not the case: only 41% (146/354) of IRAC color 
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Figure 4. 24 luminosity distribution of the star-forming and AGN sub-samples at 0.05 < z < 0.65. Left: the rest-frame L{24 ^m) 
distribution of 3,759 star-forming galaxies (solid histogram) and 288 AGNs (dashed histogram) identified by their characteristic power-law 
SED (based on power-law fitting and IRAC colors) and X-ray detections that pass all targeting and quality cuts. Right: the fraction of 
galaxies identified as AGN in our sample rises with redshift at 2 < 0.6. This is partly due to the fact that AGNs are intrinsically more 
luminous, as shown in the left panel. ULIRGs (galaxies with L{IR) > lO^'^ Lq or L(24 ^m) > f .5 X 10^^ Lq) are dominated by galaxies 
exhibiting AGN activity that are excluded from our LPs of star-forming galaxies. 

selected AGN exhibit identifiable power-law character- 
istics. However, Figure |3] (right) illustrates that there 
are a few galaxies in the power-law selected AGN locus 
but not identified as AGN by t he MIR power-law fit- 
ting. These galaxies satisfy the S tern et al. (2005) crite- 
ria, which suggests that adopting the Istern et all (l2005h 
AGN selection criteria helps augment the completeness of 
power-law AGN selection in our case. Among these 354 
IRAC color selected AGN, 159 are at 0.05 < z < 0.65 
and pass all targeting and quality cuts, and are detected 
at 24 /zm. 

As another test for MIR-emitting AGN, we used the 
behavior of the IRAC [5.8] - [3.6] color. For z > 0.2, 
both bands should be on the Rayleigh- Jeans tail of the 
photospheric SED without substantial contamination by 
infrared excess emission from normal star formation (see, 
e.g., the iRiekc ct al. (2009) SED templates). However, 
the warm dust around an AGN should appear as an 
anomalously red color. Many of the objects identified 
in this way are still dominated by aromatic emission, as 
judged by their full set of IRAC colors. However, a small 
fraction (21 objects at 0.05 < z < 0.65) may have IR 
excesses associated with AGN that are not identified by 
other methods. These 21 objects amount to only 7% of 
our AGN sub-sample. Without other independent AGN 
signatures we cannot determine whether the 24 micron 
emission of these galaxies is dominated by AGN. Our re- 
sults are independent, within the errors, of whether we 
include or exclude these sources from our star-forming 
galaxy sample. 



3.2.4. Summary of MuHiwavelength AGN Selection 

We found that optical line diagnostics (the BPT di- 
agram) select AGN that frequently have 24 /xm fluxes 
dominated by emission from star formation (Section 
I3.2.ip and therefore that these diagnostics should not 
be used as a basis to exclude AGN from the MIR 



star-forming galaxy sample. Methods based on power- 
l aw characteristic of AGN (power-law SE P selection 

J Donlev et al.l 120081 ) and color criteria, e.g. iStern et all 
2005D) yield AGN s amples more relevant to our 24 fim 
study (Section [3231) • On the [3.6]-[4.5] vs. [5.8]-[8.0] 
color-color diagram. X-ray selected AGN occupy regions 
encompassing the locus of power-law selected AGN and 
extend t oward the bluer [5.8]— [8.0] and [3.6] — [4.5] c olor 
regions (jCardamone et al.|[2008t iGoriian et al.ll2008D . in- 
dicative of older stellar populations, suggesting that X- 
ray selection helps augment the completeness of AGN 
selection by the MIR selection methods. Therefore we 
combine MIR selection of AGN with X-ray selection to 
create a list of AGN to be excluded from our 24 /im LF 
study. 

The Venn diagram (Figure [21 right) illustrate the sig- 
nificant overlapping o f AGN selec t ed by X-ray, Power- 
law fitting, and the IStern et aLl ()2005f ) criteria. At 
0.0 < z < 0.65, we find a total of 712 AGNs, 288 of 
which are detected at 24 /j,m, pass all targeting and qual- 
ity cuts, and are in the redshift range 0.05 < z < 0.65 
where we use to construct LFs. Figure [H illustrates the 
fraction of AGN as a function of 24 /im luminosity and 
redshift. Note that objects with L(24 /xm) > 10^^ Lq 
are mostly AGN. Most of the trend in the AGN fraction 
with redshift seen in Figure are a consequence of the 
increasing fraction of AGN at higher L(24 /im). Higher 
redshift bins are dominated by higher luminosity galax- 
ies because of the survey flux limit and hence will show 
higher AGN fractions. 

Excluding these 288 AGNs from the 4,047 galaxies at 
0.05 < z < 0.65 that are detected at 24 /im, we obtain 
a sub-sample of 3,759 star-forming galaxies to construct 
the LFs. 



3.3. Methodology 
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We estimate our LF at redshift 0.05 < z < 0.65 using 
the l /Vmax formahsm ()Schmidtlll968tlHuchra fc SargentI 
Il973f) . and the local 24 /xm LF (hereafter local LF or 
LLF) using the parametric maximum likelihood estima- 
tor (MLE; e.g . . ISandaee et al.l [19791: lYahil et al.l [T991 
[Marshall et~al1 119831 ). The local LF is then evolved to 
fit the higher redshift l/Vmax LFs using a minimiza- 
tion to determine the evolution of the LF as a function 
of redshift (Section 13. 5|) . 

The l/Vmax estimator is non-parametric in the sense 
that it does not assume a functional form for the LF. 
However, it is sensitive to density variations due to large 
scale structures. The large areal coverage of the ND- 
WFS helps to mitigate this problem. But at the lowest 
redshifts that we observed, even a 9 sq. deg field cov- 
ers a sufficiently small volume to show significant cosmic 
variance. The most prominent example is a ^^50% over- 
density near z ~ 0.15 (see Figure [T]). 

We construct the 1 /Vmax LF by first computing a spec- 
troscopic selection function that takes into account the 
parent catalog's incompleteness due to imaging and tar- 
geting issues. These issues include sparse sampling due to 
priority confiicts in Hectospec fiber assignment, incom- 
pleteness due to fiber allocation limits (such as the prox- 
imity limit), and redshift extraction failure (Section l2.2p . 
The selection function is then applied to each galaxy as 
a weight, Wi, that corrects for incompleteness. The Vmax 
in our analysis is taken to be the volume enclosed by 
the maximum redshift that allows an object to be in- 
cluded in both our I-band (spectroscopic targeting) and 
24 fim samples. That is, we take the smaller of the red- 
shift at which F^{2A ^m) would drop below 0.27 mJy 
or the redshift at which the object will have Ic > 20 
mag. We estimate the 1/Vmax LF in four independent 
redshift bins ranging from redshift z — 0.05 to z = 0.65 
with a bin width of 0.15. The lowest redshift region (e.g. 
z < 0.05) was excluded to avoid potentially unrealistic 
number densities due to the very small volume and num- 
ber of galaxies and also to minimize the effect of pho- 
tometric bias due to nearby, large galaxies that extend 
beyond our 24 fim PSF. This redshift cut excludes only 
82 objects (2% of the sample). In each bin, we take Vmin 
to be the volume enclosed by the Zmax of the adjacent 
lower redshift bin. 

The MLE method was used to determine the shape of 
the local 24 /im luminosity function (LLF). MLE requires 
no binning of data and is insensitive to density variations, 
which is ideal for the lowest redshift regime where large 
scale structures affect the 1/Vmax LF the most. We pa- 
ra meterize the LF w ith the double-power law proposed 
bv lYahil et all ([1991 . 



*(> L)^C 
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L 
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which corresponds to the differential luminosity function 

(2) 



L L + 



This functional form is commonly used to de scribe the 
MIR luminosity function (S hupc et al. 1998: [Xu et al.l 
199^: 'Xu"2000': I Sanders et~S.l i2003t iMagnelli et al."'2009'). 
The .Xu (2000, ) template was also used bv .Le Floc'h et al.. 



([2005D to describe MIR LF evolution to z ~ 1. It 
shoul d be noted that the turnover of the Schechter func- 
tion ([Schechteii [1976D frequent ly used for op tical LFs 
is too steep to fit the IR LF (IRieke fc Lebof skv 198^ 
lYahil et al.l[T99T[: ILe Floc'h et al.ll2005l) ! 

We determine the parameters by maximizing the likeli- 
hood function over the parameter space of the bright-end 
slope, a, the faint-end slope, /?, and the characteristic lu- 
minosity, L*, using the sample members in the redshift 
range of 0.05 < z < 0.25. The normalization factor, C, 
has to be determined separately (since the factor can- 
cels and does not appear in the likelihood function) by 
matching the integrated number density expected from 
MLE to that observed in the l/Vmax LF of the same 
sample. We constructed the LLF for both star-forming 
and all-galaxy samples. 

For the purpose of constructing a local luminosity func- 
tion, the redshift range 0.05 < z < 0.25 cannot be con- 
sidered purely local. Especially with a population with 
the strong luminosity evolution of IR galaxies, the L, es- 
timate from such a wide redshift range would inherently 
be brighter than the true local L». To illustrate this, if 
we assume that the luminosity of IR galaxies evolves as 
{1 + z)^'^ and the median redshift of the sample used for 
the LLF is z = 0.15, we would expect a 60% luminosity 
increase from the actual Li^{z = 0). There are a number 
of possible approaches to solve this problem. The first is 
to fit the luminosity evolution as a parameter in the MLE 
simultaneously while constraining the shape parameters. 
In practice such a fit would be difficult due to the addi- 
tional free parameters and the degeneracy between lumi- 
nosity and density evolution. One might need to fix some 
shape parameters (such as the faint-end slope) to achieve 
a plausible fit. The second approach would be to limit 
the redshift range of the sample to estimate the LLF. 
However, reducing the redshift bin size would compro- 
mise the constraints on the bright-end slope; we need to 
estimate the LLF using a sample selected to z < 0.25 to 
encompass enough volume to find rare, bri ght galaxies. 
Other s hav e fa c ed th e same problem. The IShupe et al.l 
(|1998[ ) and [X^ ([20001 LLFs contain few bright galaxies 
as far as z ~ 0.2. Although their median redshift is much 
closer to local, the luminosity evolution could result in 
a more luminous L* than the true local value. Another 
way to limit the redshift range to the local volume would 
be to rely on an MIR sample from all-sky IR surveys 
with shallow limiting fluxes such that the sample is lim- 
ited to very nearby galaxies. But since all-sky surveys 
lack complete spectroscopic redshift follow-up and usu- 
ally have to compile redshifts from several catalogs, the 
spectroscopic selection function of the combined sample 
may not be well understood. Robustly constraining a 
true local L, is important because the local L* is often 
used to constrain the evolution of the LF at higher red- 
shifts, such as in [Le Floc'h et al.l ([2005D . who used the 
[Sanders et all ((2003D local L* derived from the IRAS 60 
/xm-selected sample. 

Our solution to this problem is to take advantage of 
the evolution information in our full sample, an option 
not available to studies limited to z < 0.2. The evolu- 
tion constraint from the higher redshifts can be used to 
constrain the evolution and extrapolate L* back to the 
true value at z = 0.0 (Section [33[) . 
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Table 1 

Parameters of the Mid-IR Local Luminosity Function in Comparison with Other 

Works 



Reference'^ 


a 




L*(24 fim) 
(IO^Lq) 


C 

(lO-^Mpc-^) 


All Galaxies 
Star-Forming 
ShuDe et al. (1998) 
Xu et al. (19987 
Sanders et al. (■20031 


0.36 ± 0.04 
0.37 ± 0.04 
0.437 ± 0.032 
0.47 ± 0.13 
0.6 ± 0.1 


2.17 ± 0.35 
2.36 ± 0.41 
1.749 ± 0.067 
2.20 ± 0.13 
2.2 ± 0.1 


4.07 ± 0.67 
4.27 ± 0.71 
4.47 ± 0.40 
6.40 ± 2.25 
4.4 


1.2 ± 0.8 
1.2 ± 0.8 
1.0 ± 0.6 
0.9 ± 0.5 
0.8'' 



^ Results from other studies are corrected for their bandpass differences to reflect 
estimates at 24 /xm using the IRieke et all ( 2009) SED library. Redshift ranges for 
each study are the following. This study: 0.05 < z < 0.25, median redshift of 0.17 in 
both samples, with L* evolve back to 2 = 0.0 (for both al l-ga laxy and star-forming 
sample) simultane ously with the e volution fit (see Section 13.51 1: Shupc ct al. (1998}: 
0.001 < 2 < 0.07; IXu et al.l (WM) : 0.02 < z < 0.22; iSanders et al.i (12003) : median 
and mean redshifts of 0.0082 and 0.0126, respectively. 

^ ISanders et al.l II2003I) do not report their normalization explicitly but a power-law 
break fit to their 1/Vmax LLF (see their Table 6) suggests a value of 8 X 10~* Mpc~^. 
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Figure 5. 24 fim and I-band bivariate luminosity functions (BLF) 
were constructed to quantify the effects of incompleteness at 24 
ixm of MIR-bright, optically-faint galaxies (the bivariate selection 
bias). The BLF construction to quantify the bivariate selection 
bias for the sub-sample at redshift 0.35 < 2 < 0.5 is illustrated 
in this figure. The blue dots illustrate the sample observed at 
0.20 < 2 < 0.35 evolved in both I-band and 24 ^m luminosities 
to match the distribution at 0.35 < 2 < 0.5 (shown in red dots 
as a comparison). The solid vertical line is the luminosity, corre- 
sponding to the I-band limiting magnitude (/ = 20) at 2 = 0.35, 
below which a galaxy will be excluded from our 0.35 < 2 < 0.5 
sub-sample. The values indicated in the missing fraction column 
are the percent of the 24 ^m LF potentially excluded due to the bi- 
variate selection bias (fraction of 24 fj,m LF contributed by objects 
represented by blue dots on the left of the thick vertical line) in 
each L(24 fim) bin (separated by the horizontal dotted lines). The 
horizontal dashed line indicates the luminosity corresponding to 
the 24 ^m completeness limit (0.27 mjy) at 2 = 0.50, above which 
24 fira sub-sample at 0.35 < 2 < 0.5 is complete. We quantified 
the missing LF due to this effect to be less than 2% at L(24 fira) 
above the completeness limit, indicating that our spectroscopic tar- 
geting method did not exclude significant numbers of MIR-bright, 
optically-faint galaxies. 

3.4. Bivariate Selection Bias 

Since AGES spectroscopy was targeted on a flux- 
limited I-band parent catalog, our sample could be af- 
fected by biases created by excluding 24 /xm- bright, op- 
tically faint galaxies. This bivariate selection bias may 
result in a sample that systematically excludes galaxies 



with optically obscured star formation. 

We quantify the possible bias by estimating the frac- 
tion of the 24 ^m number densities excluded because 
their contributing sources fall below the I-band flux limit. 
This can be done by constructing a 24 fiin LF as a func- 
tion of I-band luminosity (bivariate luminosity function, 
or BLF) of a model sample constructed by evolving a 
sample at lower redshift, where it is complete to a fainter 
I-band luminosity, to a higher redshift bin of interest. To 
illustrate this experiment, consider the BLF in Figure [5j 
In this case, we are investigating the bivariate selection 
effect for 0.35 < z < 0.50. We construct a model sam- 
ple by evolving a sample from 0.20 < z < 0.35 (which 
is complete to deeper luminosity) in both I-band and 24 
/im to construct a BLF at 0.35 < z < 0.50. As illustrated 
in Figure [5l some objects would fall below the luminosity 
corresponding to the I-band flux limit (/ — 20.0) at the 
low-redshift edge of the bin (z = 0.35) and would there- 
fore be excluded from our sample, resulting in a fraction 
of the missing 24 fj,m LF. We found, however, that this 
fraction is small (< 2%) in this redshift bin, for 24 fim 
luminosities above our 24 /xm completeness flux. 

To model the evolu tion of L(I) w e refer to 
iWillmer et all (|2006D and lBrown et al.l (|2007| ) who found 
that the optical characteristic magnitudes brighten by 
~ 1.3 mag per unit redshift, which would correspond to a 
luminosity evolution law Lcsc (1 + z)^'^ atO < z < 1. For 
the 24 ^m luminosit y evolution i nodel , we adopted L oc 
(1 + zf -^ based on IPozzi et al.l (l200l : iLe Floc'h et"all 
(2005!) and lPerez-Gonzalez et all 120051 ). 

We repeated this test for the other two higher redshift 
bins and found that the bivariate selection bias is neg- 
ligible in the 0.20 < z < 0.35 bin; none of galaxies is 
excluded from the sample in this redshift bin due to bi- 
variate selection bias. Only in the highest redshift bin, 
0.50 < z < 0.65, does the missing fraction become more 
significant. In this bin we found the missing fraction in 
the L(24 ^m) bins of lO^O-^, 10l°■^ and 10" i Lq to be 
8.7%, 4.4%, and 0.0% respectively. 

The bivariate selection biases are under 5% for all but 
the lowest luminosity bin of our highest redshift sub- 
sample. We correct for it by by applying a correction 
factor based on the missing fraction estimated for each 
bin of the l/Vmax LF. This correction is included in the 
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Figure 6. Mid-IR local luminosity function (LLF). We construct the LLF for the all-galaxy (solid black line) and star-forming (dashed 
red line) samples. The 1/Vmax LLF points are shown in black squares above the luminosity corresponding to the 24 /xm completeness 
flux limit at the upper limit of redshift range where we select the sub-sample to construct the LLF (0.27 mjy at z = 0.25, vertical line) 
and lower-limits (upward triangles) below this completeness threshold. Note that this completeness limit does not affect the maximum 
likelihood estimator (MLE) LLF. Uncertainties in the bri ght and fa i nt-en d slopes of the MLE LLF derived from the jackknife analysis are 
illustrated by the shaded area. The mid-IR LLFs from iRush et al.l ||T993|) : [gtiupe et al.l llT998l) and [Sanders et al.i (i200a ) are corrected for 
bandpass differences and shown in circles, dot-dashed line, and stars, respectively. The consistency of our LLF shape with previous work 
shows that it is a reliable local anchor point to study the evolution of the LF. 




log[L(24 ^Lm)/l^] log[L(24 

Figure 7. Evolution of the 24 /im luminosity function (LF) of star-forming galaxies (black circles) and all galaxies (gray circles, error bars 
omitted for clarity) for four redshift bins between 0.05 < z < 0.65. The dotted LF is the local 24 ^ra LF of star-forming galaxies determined 
by the maximum likelihood estimator and the s oli d line in each panel is the local LF evolve to fit the star- forming galaxy 1/Vmax LF of 
the corresponding redshift bin (see Sections 13. 51 and 14. 211 . The vertical dot-dashed lines indicate the luminosity corresponding to the 24 /^m 
completeness flux at the upper limit of the redshift bin, above which our sample is complete. We note that the f/Vmax LF in the lowest 
redshift bin could be affected by the overdense region near z ~ 0.15 (see Figure [!}. The strong evolution of the LF can be described by 
pure-luminosity evolution. 



EVOLUTION OF THE SFR OF GALAXIES AT 0.0 < z < 1.2 



11 



LF shown in Figure [71 In a broader picture, this test in- 
dicates that our spectroscopic targeting method does not 
systematicaUy exclude MIR-bright, optically-faint galax- 
ies from the final sample. 

3.5. Evolution of the Luminosity Function 

The evolution of the 24 /Ltm LF is quantified using two 
methods. The first is to evolve the local LF template 
(Section [3^ to fit the 1/Vmax LF at higher redshifts so 
we can estimate the characteristic luminosity, L*, and 
any evolution in comoving number density as a function 
of redshift. The second is to determine the luminosity in 
each redshift bin at which the cumulative number density 
reaches a common threshold. 

In the first method, we parameterize the LF evolution 
assuming a functional form of 

iP{L,z)=giz)^[L/fiz)] (3) 

where 4>{x) is given by Equation (1). The terms f{z) = 
(1 -I- z)"^, and g{z) (1 + z)"" describe the evolution 
of the luminosity and number density with redshift, re- 
spectively. We then take our LLF determined by the 
MLE method (Section 13.31) as a local template and fit it 
to the 1/Vmax LF in each higher redshift bin using 
minimization and the median redshift of the bin. By us- 
ing the LLF as a template, we assume a fixed LF shape 
across the redshift range, and focus on the change in 
L, and C as a function of redshift. To fit the LF tem- 
plate to the 1/Vmax LF, we only consider 1/Vmax data 
points above the luminosity corresponding to the 80% 
completeness flux threshold (0.27 mJy) at the upper red- 
shift limit of the bin. Below this luminosity our sample 
in the corresponding redshift bin is no longer complete. 
The fit is done in two ways: fitting all redshift bins si- 
multaneously to estimate the global evolution of the LF 
at 0.0 < z < 0.6, and fitting each bin individually to ex- 
plore any bin-specific deviations from the global law. We 
also integrate the individually fitted LFs over luminosity 
to determine the luminosity density, which will be used 
to estimate the star formation rate density (p. Section 

EH). 

As mentioned in Section l373l the LLF is affected by the 
evolution of IR galaxies within the redshift bin where the 
LLF was constructed. With a known median redshift of 
the LLF sample (z — 0.17 for both star- forming and all- 
galaxy samples) , we can use minimization to give the 
best fit evolution law that would evolve the L* of the LLF 
from L^{zmedian) back to i*(z = 0) and simultaneously 
evolve the LLF with L*(z = 0) to fit the 1/Vmax LF at 
all higher redshift bins. After using the global luminosity 
evolution constraint to evolve L* back to z = 0, we can 
then use L*(z = 0) as a local anchor point to quantify 
the luminosity evolution in each bin individually. 

In the second method to quantify luminosity evolution, 
we determined the luminosity at which the correspond- 
ing cumulative number density reaches a common arbi- 
trary threshold in each redshift bin. The individual fit 
of the LF is integrated downward in luminosity from the 
bright-end until the integrated number density reaches 
10^'^Mpc^'^. This method robustly quantifies the lumi- 
nosity evolution because the bright-end of our LF is very 
well constrained and little affected by incompleteness. 



3.6. Estimation of SFR 

We adopt the SFR calibration of lRieke et all ()2009D to 
estimate the SFR from the observed F^{2A ^m). Every 
SFR indicator using infrared luminosity is based on the 
assumption that UV radiation from young stars is effi- 
ciently reprocessed into IR radiation, and that the frac- 
tion of UV radiation escaping the star-forming regions 
without being reprocessed is small. This escape fraction 
depends strongly on the total-IR luminosity since more 
luminous IR emission implies greater UV and optical ex- 
tinction. This relationship has been studi ed by compar- 
ing the L{UV) and L( TIR) of galaxies (jSchmitt et al.l 
120061; iBuat et al.ll2007f ). who found that the UV con- 
tribution to the total luminosity (the UV escape frac- 
tion) is about 20% for galaxies with L(TIR) equivalent to 
L{2A fim) ~ decreases rapidly to 8% and 

2.5% at i(24 ^m) ~ lO^-^ Lq and i(24 ^m) ~ lO^" ! 
Lq, respectively. Therefore the calorimetric condition 
can be assumed for our normal star-forming galaxies, 
whose L(24 /im) are in the range of 10^ — 10^^° Lq, al- 
though the error will increase between 10^ and 10*-^^ 
Lq. This is relevant only at t he lowest redshifts. We 
note that the iRieke et all (|2009[ ) 24 ^m SFR calibration 
has a scatter of < 0.2 dex, but a relative trend such as 
the evolution law will not be affected by this scatter 
so long as the sample is large. The luminosity density 
that we use to estimate the star formation rate density 
(p*) is calculated by integrating the fitted LFs (Section 
13.51) in each redshift bin over luminosity from 10^ to 10^^ 

The IR SFR indicators are calibrated with local SEDs. 
Employing these indicators at higher redshifts carries 
an assumption of no or little SED evolution over the 
relevant r ange of z (to 1 . 2). This assumption i s sup - 
ported by MarciUac et al. (20061 iBavouzet et al.l ()2008f ) 
and Magnelli et a l. (2009). There are indications that, at 
redshifts of z ~ 2, the calibration of the observed 24 pm 
fiux density and SFR r nay shift from values deduced from 
local SEP templates (iPaoovich et al.l 120071 : iRigbv et all 
120081: iFarrah et all I2008D . The primarv cause of these 
changes appears to be that the aromatic bands tend to 
be stronger for a given luminosity at high redshift, so 
the SEDs resemble more closely those of local galaxies 
with significantly smaller total infrared luminosity. The 
probable cause of this behavior is that local ULIRGs 
tend to have intense star formation concentrated into 
very optically thick nuclear regions, whereas galaxies 
of similar infrared luminosity at high redshift proba- 
bly have more widely distributed star forming regions of 
lower optical depth. It is unlikely that this effect signifi- 
cantly alters the calibration of SFRs in this paper. First, 
the effect appears to be weaker for normal star-forming 
galaxies and galaxies of LIRG- level luminosity even at 
high redshift ( Rigbv et al.ll2008l) : the galaxies that dom- 
inate our study are of lower luminosity still. Second, 
given their relatively low redshifts, the observed 24 /zm 
fiux densities from our galaxies are not strongly affected 
by the aromatic band strengths. Third, the evolution- 
ary state of our galaxies should resemble local galaxies 
at least as much as it does those at z ~ 2; in fact, 
iMagnelli et"all ()2009[) show that the SEDs of infrared 
galaxies for 0.4 < z < 1.2 resemble those of local galax- 
ies reasonably closely. 
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3.7. Uncertainties Analysis for the LF and SFR 

We use jackknife resampling to estimate the uncertain- 
ties in our l/Vmax LFs and the LLF parameters (a, (3, 
L*, and the normalization factor; see Section r3.3|) . as well 
as star formation rate density and evolution laws derived 
from the LFs. 

The AGES survey used 15 overlapping tiles corre- 
sponding to the Hectospec field of view. These tiles are a 
natural way to build 15 jackknife samples, each dropping 
one of these roughly 1-degree diameter tiles and then us- 
ing the variance in parameter estimates for these 15 sub- 
samples to estimate the uncertainties. This procedure 
will include Poisson uncertainties and sub-degree-scale 
cosmic variance effects. It will underestimate uncertain- 
ties due to larger scale correlations between galaxies. 

We can estimate the large scale effects by shifting the 
entire LF estimation for each redshift bin by our es- 
timate of the total cosmic variance and recalculating 
the parameters. The cosmic variance on the scale of 
our survey volum e (~ lO^Mpc^) can be estimated us- 
ing [Peeblei (|1980( ) and by assuming the power- law corre- 
lation function estimated from the SDSS for late-type 
galaxies, £ = (r/rn)~\ wit h ro = 3.6 h^^Mpc and 
7 = 1.7 (jZehavi et all l2005| ). In the bins of redshift 
0.05 < 2 < 0.20,0.20 < z < 0.35,0.35 < z < 0.50, 
and 0.50 < z < 0.65, we expect the uncertainties on the 
LF due to cosmi c variances to b e 15%, 11%, 10%, and 
9%, respectively (|Dai et al.ll2009t ). 

4. RESULTS 

Here we will discuss our determination of the local 24 
/im LF and compare it to other studies (Section 14.11) 
to establish a baseline for our study of LF evolution to 
z ^ 0.6 (Section 14. 2|) . Next, we discuss the evolution 
constraints based on the combination of our Boote s field 
data and the 35 /im LF evolution from Magnelli et al.l 
|2009) who focused on 0.4 < z < 1.3 (Section li^]) using 
the FIDEL survey. The combination of these two studies 
lets us determine the evolution from < z < 1.2 in a 
unified and accurate way. We will then present the local 
star formation rate density and the star formation rate 
density evolution derived from the LFs fSection l4.4|) . 

4.1. Local 24 fJ-m Luminosity Functions 

The local 24 /im LFs for both the star-forming and all- 
galaxy populations are shown in Figure [51 Both LLFs are 
constructed from a sample at 0.05 < z < 0.25 (median 
redshift z = 0.17) and evolved back to z = 0.0 using 
the global evolution law estimated in Section 13.51 Their 
shape parameters are tabulated in Table [T] along with 
other estimates of the LLF using the double-power law 
form of equation. 

The star-forming and all-galaxy LLFs have identical 
faint-end slopes, but the star-forming LLF has a slightly 
steeper bright-end slope and a brighter characteristic lu- 
minosity. Both differences are within the uncertainties of 
the respective parameters, but nevertheless the steeper 
slope of the star-forming LLF is expected because the 
sample excluded the IR-luminous AGN, as illustrated in 
Figure [H This comparison shows that the overall im- 
pact of AGN on the LLF is small, and that any residual 
unidentified AGN should have little effect on our conclu- 
sions. In Figure |6] we present our LLF along with the 



Table 2 

The 24 fim and Total-IR Luminosities 
Above Which the Cumulative Number 
Density of Star-Forming Galaxies is 
10-*Mpc-3 for Each Redshift Bin 



Redshift Bin 


L(24 ^jm) 


L(TIR)=' 




IOIOL0 


IOIOL0 


0.05 < 2 < 0.20 


1.2 ± 0.4 


9.2 ± 2.9 


0.20 < z < 0.35 


1.7 ± 0.4 


13.2 ± 3.0 


0.35 < 2 < 0.50 


2.5 ± 0.5 


18.9 ± 3.7 


0.50 < 2 < 0.65 


3.9 ± 0.7 


28.1 ± 5.3 


L(TIR) is 


estimated 


using the 



IRieke et al I l|2009l ') formalism. 



bandpass a nd cosi nologically adj usted 1/Vr r,.„x LLF of 
iRush et al.l (|1993D . lSanders et all (I2003D a n dlPozzi et al.l 
( j2004D . and the MLE LLF of IShuoe et al.l fOWW . 

One of the most important parameters of the LLF is 
the value of L*(z — 0). Our L» estimates for the star- 
forming and all-galaxy populations are (4.27±0.71) x 10^ 
and (4.07±0.67) x 10^ L©, respectively. We compare our 
results with other works that adopt the double power- 
law form of LF by correcting for cosmological assump- 
tions and b andpasses using the ,Rieke et al. (2009) SED 
library. The IShupe eFal] (|l99l LLF using IRAS 25 fj.m 
measurements reports a local L* of (4.67 ± 0.41) x 10^ 
Lq, which corresponds to (4.47 ± 0.40) x 10^ Lq when 
corrected for the differenc e between the IRAS 25 /im and 
MIPS 24 fim bandpasses. iXu et al.l 1(1998,) report an ISO 
15 /im local L* of (4.79 ± 1.69) x 10^ Lq, corresponding 
to (6.40 ± 2.25) X 10^ L© at 24 /im. ISanders et all (pOOl 
report a local total-IR L* of IO^^'^Lq, which corresponds 
to a 24 /im L* of 4.4 x 10^ Lq. Our estimate of L* agrees 
well with all these other studies. More importantly, the 
extrapolation of the LLF L^, back to z = 0.0 using the 
luminosity evolution constraint from our full sample en- 
sures that these values reflect an accurate value at z = 
in addition to providing the best-match LLF template to 
study the evolution of the LF in our sample. 

4.2. Evolution of the 24 ^im LF 

Strong luminosity evolution is evident in the 1/Vmax 
LFs over the redshift range 0.05 < z < 0.65 shown in Fig- 
ure [T] We considered two evolutionary scenarios for the 
24 /im LF parameterized by Equation (3), including pure 
luminosity evolution and pure density evolution. These 
scenarios were tested by evolving our local 24 /im LF to 
fit the 1/Knax LF at higher redshift by minimization. 
We cannot constrain a combined luminosity and density 
evolution model because we do not sample below L* in 
these higher redshifts sub-samples. The density evolu- 
tion model can be rejected on the basis of the values 
with high confidence. Other studies have fou nd that the 
density evolution is small. ILe Floc'h et al.l ([2005 ) and 
IPerez-Gonzalez et al.l (j2005() derived a density evolution 

laws of oc (l-Hz)°- ^-°.6 smAcj)* (x (l-Hz)^ -°^°-^ respec- 
tively. More recently, IMagnelli et al.l ()2009l ) used a stack- 
ing analysis to study the evolution of the 35 /zm LF and 
reported density evolution given by (p* (x {1 + 2;)O i±o 'i' 
at z < 1.2. This implies a density variation of about 
5% ± 40% at 0.0 < z < 0.6, which is comparable to the 
uncertainties due to cosmic variance in the density. 
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Redshift (z) 



Figure 8. Evolution of the characteristic 24 /im luminosity of galaxies, L*, to 2 ~ 1.2 based on our results (red circles) and that of 
IMagnelli et a l. (2009) using FIDEL (red squares). The light shaded area is th e luminosity e volu tion law L^, oc (1 + z) "^ bounded by ai^ 
of 2.7 and 4.2 from the lowest and highest uncertainty limits of iPerez-Gonzalez et al.l 1)20051 ^ and IMagnelli et al.l II2009I ). respectively. The 
combined evolution law, combined oc (1 + z)^ *^'^ '^, is shown as a solid line with the uncertainties from the combined fit shown by the 
dark shaded area. The local point for the global fit (red star) is that of IShupe et al.l II1998I ) converted to 24 ^ra. Two data points were 
excluded from the fit (see details in Section |4.3| I and are shown in gray. Inset: the combined of luminosity evolution (solid line) with the 
same symbols as the main figure (error bars arc omitted for clarity) and the fit (dotted line) based on the Bootes data alone, (1 + z)^'*^"'^ . 




0.2 0.3 0.4 

Redshift log(l +z) 



Figure 9. Evolution of the star formation rate density (p») as a functi on of redshift assuming lSalpeterl | |1953 '| IMF. Our p* estimates are 
shown in large black squares and our extrapolated local point (Section 14.111 shown in gray, compared to the p* from the following. The 
extinction-corrected UV p, from the compilation of Hopkins ( 2004) and Hopkins & Beacoinl H200fl ) are shown in blue circles; IR p» in red 
squares (Hopkins 2004; Perez-Gonzalez et al. 2005; Reddv et al. 2008); gray dots ar e the radio ( 1.4 GHz), X-ray, and extinction-corrected 
[OH], Ho, and H,3 p, (Hopkins 2004; Seymour et al. 2008; Smolcic et al. 2009; Zhu et al.l[2009t) . The agreement of our p* normalization 
and those from extinction-corrected optical and UV observations confirms that a majority of star formation at z < 1.2 is obscured and can 
be directly traced by 24 pm observation. 
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Table 3 

Evolution of the Characteristic 24 /im Luminosity, 
L«{24 ^m), and the Star Formation Rate Density, p* 



to z 


~ 1.2 




Redshift 


L*(24 pm) 


log{p.)^ 




109 Lq 




0.00 (this work) 


4.27 ± 0.71 


-1.83± 0.08 


0.00 fShuDe et al. 19981 


4.47 ± 0.40 


-1.71± 0.04 


0.05 < 2 < 0.20 


7.25 ± 1.12 


-1.60± 0.07 


0.20 < z< 0.35 


10.6 ± 1.8 


-1.43± 0.08 


0.35 < z < 0.50 


15.5 ± 2.4 


~1.26± 0.08 


0.50 < 2 < 0.65 


23.6 ± 3.6 


-1.07± 0.07 


0.7 < 2 < l.O*" 


32.8 ± 7.8 


-0.92± 0.11 


1.0 < 2 < 1.3*" 


65.8 ± 18.8 


-0.61± 0.14 



p, estimates assume a lSalpet"e3 1)19551 ) IMF and are 
given in the unit of M0yr~^Mpc~^. Assuming the 
IKroupal SWO^ ) IMF would yield 0.66 p* given here. 

^ Based o n ^FIDEL luminosity functions 

UMagnelli et al.ll2009l '). 

Using a pure luminosity evolution model, we obtain an 
evolution law given by L,(24 /im) oc {1 + z)^-^^'^-^ for the 
star-forming sample. This luminosity evolution is shown 
as a fit to the l/Knax LF in Figure [71 One of the evo- 
lution descriptions of iPerez-Gonzalez et al.l (120051) uses 
pure luminosity evolution with the lSchechteii ()1976[ ) form 
of LF. They report luminosity evolution of (1 -I- 2)3.6±o.3 
for z < 0.8, agreeing with our result. The evolution of 
the LF was also quantified by determining the luminos- 
ity above which the cumulative number density reached 
IQ-^Mpc"^; as tabulated in Table H This luminosity 
increases by more than a factor of three from a redshift 
of 0.05 to 0.65. We should note that while we cannot 
jointly constrain luminosity and density evolution, the 
simple power-law model is an adequate description of 
our data. Improving on it would require a deeper sample 
in a comparably large area to the Bootes field. 

Our method to determine evolution of the LF assumes 
no evolution in LF shape from the local template. The 
template, however, yields a better fit at low redshift than 
in the higher redshift bins. Discrepancies between the 
template and the data are most visible at 0.50 < z < 
0.65 (Figure [T])- This discrepancy could be due to an 
intrinsic evolution toward a shallower bright-end slope at 
higher redshifts. On the other hand, it could indicate a 
contribution from unidentified AGN in the higher 24 fim 
luminosity regime. The latter explanation is consistent 
with the observation that AGN are progressively more 
abundant at higher redshifts (Figure |4]) as well as from 
Figure [7] itself, which shows that the bright-end slope of 
the l/Vmax LF of all galaxies (including AGN) becomes 
shallower to a greater extent compared to that of the 
star-forming galaxy LF. 

4.3. Extending the LF evolution constraint to z ^ 1.2 

We augment our LF evolution constra int at low and 
interm ediate redshifts with the result of Ma gnelli et alj 
(I2OO9I) at 35 /ir n. Th is is a useful combination because 
iMagnelli et al.l (|2009l ) used the deep (50 fiJy and 3 mJy 
at 24 and 70 /im, respectively), narrow-field (0.5 sq. 
deg.), FIDEL survey. The FIDEL sample has - 7000 
and ^ 600 star-forming galaxies detected at 24 /im and 
70 /im, respectively, with known redshifts. AGN were 
excluded on the basis of deep (200 ks - 2 Ms) X-ray ob- 



servations. Additionally, IMagnelli et al.l ()2009[ ) employ 
a stacking analysis to provide further constraints for the 
LFs. Their approach provides one of the best evolution 
constraints for the LFs at higher redshifts to date but 
will be affected by cosmic variance at intermediate red- 
shifts, where our results should be more reliable. In ad- 
dition, the rest frame 35 /im values can be corrected to 
24 /im with minimal uncertainties due to luminosity con- 
versions. 

We need to correct IMagnelli et al.l ()2009[ ) for (1) band- 
pass differences; and (2) k-correction differences (SED 
assumpti ons). To co r rect fo r bandpass differences, we 
use the iRieke et all ()2009D SED library to calculate 
a luminosity-d e pende nt conversion factor between the 
iMagnelh etalj (|2009D rest-frame 35 /im, which was ob- 
served at 70 /tm and k-corrected to z = 1 (the ob- 
served 70 /im band probes the 35 /im rest-frame SED 
at z = 1). This conversion factor varies slowly from 0.5 
at L(TIR) - IO^^Lq to 0.3 at L( TIR) ^ 10^^-^(^ - Next 
we compare t he consistency of the | Charv &: Elb azf (20011) 
SED u sed bv IMagnelh eTall (pOOa T and our[Rickc et al] 
()2009[) SED. The best test for differences lies in the 
L,(z = ) anchor poin t, L»(z = 0,35 /tm) = 10^-^^ 
Lq, th at IMagnelli et al . (2009) obtained by convert- 
ing the iShupe et al.TilTQQSil IRAS 25 /im L* to 35 /im 
(MIPS 70 /im o bserve d bandpass at z = 1) using the 
|CliariL& Elbai (|200l SED. Converting this U{z = 
0,35 /im) back to our MIPS 24 /im rest-frame (which is 
ab out 5% fainter tha n the original value at 25 /im) using 
the lRieke et"all ([2009i) SED would yield L^{z = 0, 24 /im) 
of 3. 72 X 10^ h-fL^, w hich is 20% lower than the orig- 
inal IShupe et aH (|1998D value tabulated in Table [H In 
other words, a 20% discrepancy resulted from subject- 
ing a luminosity through this conversion exercise and we 
must apply a -1-20% correction factor to the FIDEL data 
points to correct for SED differences. 

For the combined fit, we omitted two data points. 
First, we dropped the lowest redshift Magnelli et al] 
()2009t ) bin (0.4 < z < 0.7), which overlaps with our high- 
est redshift bin, because this redshift bin will be more 
strongly affected by cosmic variance since our highest 
redshift bin encompasses about an order of magnitude 
more volume. Second, we dropped our local point be- 
cause it was estimated by evolving the sample in the 
lowest redshift bin (0.05 < z < 0.20) back to z = 0.0 
based on the higher redshift bins (Section l3.5p . so includ- 
ing this p oint would "double count" the data. Instead 
we use the lShupe et al.l (|1998[ ) local point converted to 24 
/im. Note, ho wever, that our local L* estimate is consis- 
tent with Shu pe et al.l l)1998l ). The combined fit (shown 
in Figure [8]) yields a slightly shallower evolution law than 
that estimated from the Bootes field alone fSection l4.2p . 

with L^,comMned OC (l-f z)3-'^='=°-^ . 



4.4. The Evolution of the Star Formation Rate Density 

Now that we have constructed the LF and determined 
the luminosity evolution from z ~ 0.0 to 1.2 we can pro- 
vide a complete picture of the /»* evoluti on from z = 0.0 
to z - 1.2. We compute SFRs based on IKroupal (|2002[ ) 
IMF used bv IRieke et alj (|2009l ). This leads to estimates 
of a factor of 0.66 lower than those assuming a lSalpeteil 
(|1955f) IMF with a single power-law slope of —1.35 from 
0.1 to 100 Mq. 
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Our estiniates for in each redshift bin to z ~ 1.2, 
including both our re sults and the adjusted values from 
iMagneUi et all (|2009l) . are tabulated in Table |3l and 
shown in Figure |9] along with measurements from 
other wavelengths from the UV to radio. Based on the 
fit of the combined data out to z ~ 1.2, our evolution 
law is given by 

log(/9*) = (3.5 ± 0.2)log(l + z) - 1.96 ± 0.08 (4) 

At z = our fit yield s a local of (1.09 ± 0.21) x 
lO-^Mgyr-iMpc"^ for IKroupal (l200l) . which corre- 
s ponds of ( 1.65 ± 0.32) x lO-^Meyr-iMpc"^ for 

a iSalpeteii ()1955() IMF. Our estimate agrees with pre- 
vious works on the local that have properly ac- 
counted for the internal extinction of the galaxy. With- 
out such correctio ns the /j» wi l l be underestimated by 
50 - 70 % (Gallcg o et al.l[T99l PiVesse fc MaddoxHTool 



iHanish et al.. ,2006t iSalim et al.l l2007f ). The agreement 
between estimates corrected for internal extinction 
and our derived from 24 pm luminosity suggests that 
^half of the local star formation is obscured and the 
MIR provides a direct and reliable means for SFR mea- 
surement. 

Figure |9] illustrates the good agreement of the nor- 
malization of our p* estimates and those derived from 
the UV observation corrected for the effects of intrin- 
sic extinction of the target galaxies (normalization cor- 
rections can be as large as 0.7-dex; e.g.. Figure 5 of 
iSchiminovich et all (|2005l ) illustrates the exte nt of the 
requir ed correction for the UV observations). iHopkind 
(|2004| ) reports the fit to a compilation of the p* that 
has been corrected for extinction of the target galax- 
ies with an assumption of a luminosity-dependent ob- 
scuration. The resulting best fit is given by log(£*_^ = 
(3.29± 0.26)log(l-Kz)-1.98±0.04 (assuming thelKroucal 
I200I IMF), which also agrees well with our evolution 
determined from the combination of the NDWFS and the 
FIDEL results. Again, this agreement confirms that star 
formation at redshift 0.0 < z < 1.2 occurs in obscured 
environments and can be traced directly with 24 /xm ob- 
servations given a proper classification and exclusion of 
AGNs. It shoul d be noted that our re sult is consistent 
with findings of iTakeuchi et all (|2005[ ). which show an 
evolution of the fraction of obscured star formation from 
50% locally to > 80% at z ~ 1. 

5. CONCLUSIONS 

We study the evolution of 24 /im-selected galaxies by 
constructing their 24 pm (rest-frame) luminosity func- 
tions (LF). Our sample contains 4,047 galaxies with spec- 
troscopic redshifts at 0.0 < z < 0.65 in the Bootes field 
of the NOAO Deep- Wide Field Survey (NDWFS). The 
24 pm data and spectroscopic redshifts were obtained 
with Spitzer and with the AGN and Galaxies Evolution 
Survey (AGES), respectively. Our sample is a unique 
combination of wide field, spectroscopic redshifts, and 
24 pm imaging that is ideal for intermediate redshifts 
where MIR galaxy evolution has not been well-studied. 
The large area (9 sq. deg.) of the Bootes field helps 
mitigate cosmic variance, which otherwise can be a se- 
rious issue at low and intermediate redshifts. Theoreti- 
cally, the cosmic variance in our study is less than 10% 
at z > 0.2. Spectroscopic redshifts reduce cross-talk be- 
tween redshift bins and faint/bright end slope biases and 



provide an accurate number density contribution for each 
object. 

Galaxies exhibiting AGN activities characterized by X- 
ray emission and mid-IR power-law are excluded from 
our star-forming sample. We find that the optical line 
diagnostics (the BPT method) are not suitable for identi- 
fying AGN dominant at mid-IR wavelengths because the 
mid-IR emission from optically-selected AGNs is often 
dominated by star formation. We identified 288 objects 
with 24 pm emission likely to be dominated by AGN. 
Excluding them leaves a sample of 3,759 star-forming 
galaxies. 

We derive the 24 pm luminosity using calibrations 
derived from the most r ecent SED library of the mid- 
IR spectra from Spitzer (jRieke et al.ll2009( ). Our sam- 
ple is comprised mainly of normal star-forming galaxies 
(66% have i(TIR) < 10" L©) and LIRGs (33% have 
10" Lq< L(TIR) < 10^2 j^^^y Only 1% are ULIRGs 
(L(TIR) > 10^2 and we found that most of the 

ULIRGs contain AGN. 

We construct the local 24 pm luminosity function 
(LLF) as a template to study evolution of the LF. Our 
LLF was constructed at 0.05 < z < 0.25 and evolved 
back to z = simultaneously with the fit to constrain 
the global evolution of our sample. The LLF tem- 
plate for our star-forming sample is given by a double 
power-law with faint-end slope, a = 0.37 ± 0.04, bright- 
end slope, P = 2.36 ± 0.41, characteristic luminosity, 
L^z = 0) = (4.27 ± 0.71) X lO^L©, and density nor- 
malization, C = (1.2 ± 0.8) X lO^^Mpc^^. Our major 
results are as follows. 

1. The evolution of the LF at z < 0.65 can be repre- 
sented by pure luminosity evolution with the character- 
istic 24 pm luminosity of star-forming galaxies evolving 
as L*(24 pm) oc (1 -I- 2;)3-8±o.3^ demonstrate by the 
construction of the bivariate I-band and 24 pm luminos- 
ity function that the fraction of missing 24 pm-luminous, 
optically-faint galaxies due to our I-band magnitude limit 
for spectroscopic targeting is very small. 

2. We extend the constraint on the evolution of L* to 
z ~ 1.2 by c ombining our re s ults w ith the higher redshift 
results from lMagnelli et all (|2009il based on the FIDEL 
survey. The combined sample gives a slightly shallower 
luminosity evolution of L^(24 pm) cc (1 -I- z)'^ '*^°-^. The 
shallower evolution law derived by including the FIDEL 
sample at redshifts of z = 0.85 and 1.15 suggests that 
the evolution in p* may begin to slow compared to lower 
redshifts. 

3. The local star formation rate density based on our 
24 pm data i s p^ of ( 1 09 ± 0-21) x IQ-^Mgyr-iMpc"^ 
assuming the IKroupal (|2002f ) IMF, which corresponds to 
p^ of (1.65 ±0.32) X lO-2M0yr-iMpc-3 for a ISalpeted 
(|1955f ) IMF. The combined evolution constraint for the 
LF at z < 1.2 indicates that p* evolves as p* oc (1 -|- 
^■j3.5±o.2 Our normalization and evolution measure- 
ments agree well with the z < 1 extinction-corrected 
findings from other studies, which confirms that most 
of star formation is obscured and that 24 pm serves as a 
direct and reliable p* indicator in these cases. 
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